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§ . Abstract 
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' The procedure of maximization of information entropy can be used 

to improve our knowledge of parton distributions. This method has 
been applied in order to achieve improved description of the nuclear 
effect in T production due to gluon distribution modification in nuclei. 
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1 Introduction 



Information Theory was originally developed by Hartley JIJ, Nyquist 0, and 
Shannon || in order to establish a mathematical description of telecom- 
munications and to understand how information may be lost upon trans- 
mission over noisy channels. Shannon, in particular, developed a complete 
formalism in which the concept of information is quantified and important 
theorems regarding its transmission are proven. The fundamental quantity 
that measures information is information entropy. Shannon has shown that 
the information entropy is the most suitable function of the probabilities for 
emission of signals by an ergodic source that measures the magnitude of the 
receiver's uncertainty on those signals. In this sense information entropy is 
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a true measure of one's ignorance of their correct content. The larger the 
entropy the greater the uncertainty and, consequently, the information con- 
tent. Therefore, maximizing the information entropy can lead to evaluation 
of the signal probability distribution under the constraints imposed by the 
telecommunications problem at hand. 

Application of this theory in physics could be versatile. All quantum 
phenomena, for example, are stochastic in nature and are described in terms 
of probability amplitudes. Then by appropriately defining the information 
entropy of the physical system under consideration and maximizing it under 
constraints imposed by theoretical assumptions or experimental data one 
can obtain the probability amplitudes that are most consistent with one's 
ignorance of the system. This method has been used by Plastino to 
evaluate wave functions for various physical systems. It is a very powerful 
technique since it does not rely on any specific modeling but only on what 
is actually known about the system to derive best estimates for what is 
unknown. 

In this article we apply Information Theory to improve parton, specifically 
gluon, distributions in nuclei. These are suitable subjects to the method be- 
cause they are probabilistic in character. The problem whose solution we wish 
to improve is that of quarkonium suppression in proton-nucleus collisions at 
very high energies. The production of charmonium states, most notably of 
the J/ty boson as well as of bottomonium ones, especially the T resonances, 
has been observed in various experiments involving heavy nuclear targets to 
be lower that in hydrogen if the latter is multiplied by the mass number 
of the nucleus. At energies achieved at Fermilab J/\l/ and T suppression is 
very pronounced and exhibits a characteristic dependence on the momentum 
fraction of the target nucleon carried by the struck parton || . Many models 
have been developed to explain this behavior. For the purposes of demon- 
strating the Information Theory method we consider a model that is based 
on the assumption that quarks in nuclei have a finite probability to conglom- 
erate forming multi-quark color singlet states, usually called (multi) quark 
clusters |6[]. The parton distributions in such clusters differ from those in 
single nucleons and generally are concentrated to lower momentum fractions 
of the partons as the cluster becomes larger. This model supplemented by fi- 
nal state dissociation of the produced quarkonium has successfully described 
J/\I/ suppression in hadron- nucleus collisions as well as most heavy ion col- 
lisions, the latter being the topic of heated debate due to their relevance 
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with Quark-Gluon Plasma production. In this model very simple parton 
distributions have been used based on very general assumptions. However, 
the gluon distributions which play a crucial role in quarkonium production 
are poorly known in nuclei. It turns out that within this model the gluon 
distributions that solve the problem of J/\l/ suppression are inadequate to 
describe T suppression from the same experiment. We shall use Information 
Theory to improve them in a manner that maintains their applicability to 
the J/ty data and, at the same time, to enhance agreement with the T data. 



2 Information Theory 

Suppose an ergodic source of information, which can be anything from a tele- 
graphic device to a quantum system, produces signals x from some available 
ensemble X with probability distribution p(x). We define the information 
entropy of the source as || 

S = - ^2p(x, otj) In p(x, Oi), (1) 

x 

where the summation includes all instances of x in X and can indicate an 
integral if x is a continuous variable and «j are a group of fixed parameters 
in the function p. It is required that 

X>foai) = l- (2) 

X 

If the logarithm is binary then S is expressed in bits. Considered as a func- 
tional of the probability distribution, S is maximal when p is uniform; we 
know the most about a system when the signals it produces have no variabil- 
ity. To obtain the optimal function p(x, cti), i.e., to estimate the "best" set 
of parameters an we impose the extremization condition 

dS 

for all i. To ensure a maximum second order derivatives must be looked at as 
well. In this work we assume a certain functional form for p and simply want 
to determine its parameters. The procedure can be generalized to include 
cases in which we do not know the exact form of p but there are constraints 
based on data. The method of Lagrange's multipliers can be used to evaluate 
distribution functions when a number of expectation values are known [|J. 
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3 Quarkonium Suppression in p-A Collisions 



3.1 Facts and models 

The observed suppression of the J/^f, production cross section per unit 
mass number, A, in high energy hadron-nucleus j5| |7| and nucleus-nucleus |§ 
collisions exhibits a strong nuclear dependence. The hadron-nucleus data 
also show that the depletion increases with the longitudinal momentum of 
the quarkonium. These results have generated many theoretical studies. A 
variety of effects are thought to contribute and numerous models have been 
suggested. These contributions may be grouped into six major categories: 
(1) Quarkonium scattering off parton and/or hadron co- movers [|TI], [TTj]; (2) 
Glauber inelastic scattering on the nucleons fll2| ; (3) Shadowing and EMC 
distortions of the nuclear parton distributions [O]; (4) Parton scattering 



before [14| and after JT5[ the hard process; (5) Parton energy loss in the 



initial and/or the final state |16| and, (6) Intrinsic charm in nucleons [11 
Color transparency may alter the charmonium Glauber absorption JT7|. In 
the case of heavy ion collisions contribution of an unconfined state (Quark 
Gluon Plasma) has probably been very small in earlier experiments |8| but 
is currently debated due to new data []18| . In our view more than one con- 



tribution must be carefully balanced to explain the quarkonium suppression. 
Here the modifications of the initial state parton distributions (for all values 
of the Bjorken variable, x) and the final state inelastic scattering (absorption) 
are considered. 

In this article the term "EMC effect", named after the European Muon 
Collaboration, will signify any deviation from unity of the structure function 
ratio of a bound nucleon to that of a free one at any value of the variable x, 
the latter defined as the fraction of the nucleon momentum carried by the 
interacting parton. 



3.2 Quark clusters in nuclei 

The EMC effect has been studied in the framework of the expansion of a 
nuclear state \A) on a complete basis of color singlet states labeled by the 
number of (3, 6, 9 or more) valence quarks they contain. Such states, also 
referred to as multi-quark clusters, are formed when nucleons bound in a 
nucleus overlap sufficiently so that they share their constituent partons. The 
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probabilities for multi-quark cluster formation can be estimated using nu- 



clear wave functions |J or can be found by fitting DIS data JT9], P0| . The 
agreement with EMC and NMC jnj] data is excellent down to Q 2 « 2 GeV 2 
and the fit strongly constrains the quark momentum distributions in clus- 
ters. (Momentum conservation in clusters fixes the total momentum fraction 
carried by the gluons.) Good description of the small enhancement above 
unity of the EMC ratio (antishadowing) around x — 0.1 requires inclusion of 
up to 12q clusters but the essential features of the data are accommodated 
by a truncation to the 6q term and use of an effective 6q cluster probability, 
/. The observed shadowing of the structure function, F 2 , in the nucleus for 
low xb and the depletion for intermediate x combined with the QCD sum 
rules require the presence of antishadowing of Fi for other values of x. This 
antishadowing, however, need not be restricted to the range 0.05 < x < 0.2. 
In fact, this model predicts antishadowing for x > 0.8 in agreement with the 



data of Ref. [^1|] on the slope of the Ca over deuterium structure function 
ratio. In addition, the electron DIS data from SLAC confirm the excel- 
lent agreement of this model with the observed nuclear dependence for all 



x > 0.1 [22]. The model has been successfully applied to explain nuclear 



effects in Drell-Yan processes [231 and gives interesting predictions for these 
effects at RHIC energies [p4}| . Direct photon production in hadron- nucleus 
collisions |25[ has also been predicted to be altered by nuclear effects. 

The probability / depends on A approximately logarithmically and is 
0.040-0.052 for deuterium [p5| . Very dense nuclei ( 4 He) have / values larger 
than the logarithmic prediction. In the scaling limit the Nq proton-like clus- 
ter parton momentum distributions are assumed to have standard forms, 

V$\x) = B^{l-xf»\ b d N = b u N + l, (4) 
S N {x) = A N (l-x) aN , A N = x s N (l + a N ), (5) 
G N {x) = C N (l-x) CN , C N = x%(l + c N ). (6) 

The exponents that best describe the NMC data are (63 , 03, 6g, a^) = (3, 9, 10, 
11). For the gluon exponents the direct 7 data suggest (c 3 , c 6 ) = (6, 10) P5| ]. 
The valence exponents approximately follow the dimensional counting rules. 
The ocean consists of three species of quarks (and their antiparticles) with 
the strange distribution being half as large as the up (or down) quark sea 
distribution. The gluon momentum fraction is taken 5 times larger than 
the ocean one. Isospin invariance relations connect the distributions that 
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belong to the same isospin states reducing, thus, the number of independent 
paramenters. Kinematics forces the fraction of the cluster momentum carried 
by a parton in 6q clusters, x^\ to be half as large as that in nucleons, x^ 3 \ 
In this model the QCD sum rules are explicitly obeyed. Furhter description 



can be found in Refs. [G4L 25 



3.3 Theoretical cross sections 

Including 3q and 6q clusters in the nucleus the doubly differential hadron 
level cross section for quarkonium production is calculated from the a 2 order 
parton level ones (with both quark annihilation (qq) and gluon fusion (gg) 
included). The latter are convoluted with the appropriate sums of prod- 
ucts of parton distributions Hqq{xi,x\ f) and H^(xi,x^) (i = 3,6), respec- 



tively [26]. The indices 1 and 2 refer to the probe (moving in the +z direction 
in the laboratory) and the target, respectively. The duality (in effect color 
evaporation) hypothesis is applied to integrate the differential cross section 
over m 2 with m ranging from twice the c(b) quark mass, 2m c ^) > to the open 
charm (bottom) threshold, 2m,Dm)- Here = 1.864 GeV and mg = 5.278 
GeV. The duality constant, F^, is simply the portion of the total charmonium 
cross section (up to a given order) that corresponds to the quarkonium; it 
cancels in cross section ratios. The resulting cross section can be expressed as 
a function of the longitudinal momentum fraction carried by the produced 
quark- ant iquark pair, xp = 2pi,/y/s. The higher order corrections are as- 
sumed to result in a multiplicative factor, K. The transverse momentum 
dependence due to higher order diagrams is, thus, integrated over and that 
due to the intrinsic transverse momentum of the partons is neglected. These 
considerations lead to the order a 2 equation 

J—(A) Mm 2 

-KF d ° (B) dm 2 £ J® [H® a q ,(m 2 ) + H% a 99 (m 2 )] , (7) 



dxi 



/4m ,,. ■ q a 

c(b) 1=6,0 



where a q q and a gg are the partonic-level cross sections for the two hard pro- 
cesses and are the Jacobians that transform x\ and x% to xf and m 2 . 
In this way the entire quarkonium yield is found. The various final states 
being unitary rearrangements of one another have up to this point the same 
dependence on EMC-type nuclear effects at the same yfs. The perturbative 
Q 2 evolution of the parton distributions is omitted as it largely cancels in 
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the calculation of cross section ratios. In addition, the Q 2 values that are 
relevant to the production of quarkonium states are inside the scaling region 
in which the parton distributions used in this work are valid. 



3.4 Final state absorption 

After its production the cc (bb) system propagates in the nuclear medium 
developing into a quarkonium state which then decays into the observed 
lepton pairs. During this stage the system may be inelastically scattered by 
3q and 6q clusters. It is reasonable to assume that the scattering occurs with 
clusters bound in the nucleus since the break up time of the latter exceeds 
the time the pair needs to traverse the nuclear radius. The J absorption 
cross section has been measured to be a^ s ~ 3.5 mbarn/nucleon |27). This 
number has been extracted in a kinematic range in which parton distribution 
modifications are negligibly small (x% ~ 0.22, the point at which the EMC 
ratio crosses the unit line) and since it is the average absorption cross section 
over the path traversed in the nucleus, color transparency effects are already 
in it. The is attenuated a little more than the J/^f due to its larger 
radius. In this It is also assumed that the cross section for bottomonium (T 
states) absorption is of the same order as that of charmonium but smaller 
(« 3 mbarn/nucleon) due to its more compact size. 

The cross section on a 6q cluster is 2 3 / 2 times larger than that on a nucleon 
(bag model estimate) but the density of scattering centers in the medium is 
reduced when / ^ 0. Then 

(pa abs ) = v®p A [{l - f) + 2 3 / 2 /]/(l + /), (8) 

where pA is the number density of the nucleus A taken as constant within 
the nuclear volume. The average path length the cc (bb) pair travels in an 
approximately spherical nucleus of radius ta estimated by means of a simple 
geometrical (not eikonal) calculation is La ~ ItaI^- The experimentally 



measured nuclear RMS radii [£8] are employed to compute pa and L A . The 



cross section in Eq.([7|) is then attenuated by the factor 

P A = exp[-(pa abs )L A ]. (9) 
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4 Charmonium Suppression 



The nuclear dependence is extracted by taking the ratio, Ra, of the J/ty 
production cross section per unit A in collisions of a hadron, proton in this 
case, with a heavy nucleus to that in collisions with a light one and exam- 
ining its xf dependence. Since the K- factor may also depend on xf at first 
we examine the ratio of the experimental to the order ot 2 s theoretical cross 
section. In Fig. 1 we present this ratio for Cu and Be using a representative 
set of parameters for our model and the data of Ref. 0. F d is the same in 
both cases and cancels in the ratio. Clearly the i^-factor exhibits a strong 
xp dependence but, most importantly for our purposes, is the same for the 
two nuclei. This implies that it will cancel in ratios of two theoretical or 
experimental cross sections. The origin of this factor must, thus, lie in pro- 
cesses that do not depend on the size and intrinsic properties of the nucleus 
(the absorption part is included in the results of Fig. 1). The xf dependence 
of K in Fig. 1 agrees with that of the ratio of the "diffractive" to hard cross 



sections in Ref. p9| . Issues related to the relative suppression of various 
charmonium states are outside the scope of this article. 

The results for Ra are confronted with the data of Ref. || in Fig. 2. The 
dotted lines in Fig. 2 represent the prediction of the full model including 6q 
clusters and final state absorption with the lowest (highest) value of / for 
the heavy (light) nucleus and the largest values of the 6q ocean and gluon 
exponents, (ae,ce) = (12, 11); the solid ones correspond to the opposite 
/ combination and (ae,ce) = (10, 9). In order to make the influence of 
the initial and final state contributions to Ra clear in Fig. 2(d) the results 
without final state absorption are shown (short and long dash curves with 
the same connotation as the dotted and solid ones, respectively) as well as 
the result with only final state absorption, / = (dot-dash line). It is evident 
that the predictions of the full model are in agreement with the data. We 
note in passing that for large negative xp this model predicts antishadowing 
of production in p-A collisions because in this region large 22 values for 
the gluon distribution ratio are accessed. 
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5 Bottomonium Suppression 



Using the model we described earlier we can also calculate the suppression 
ratio for the T states and compare the results with the data of Ref. [f|. 
Specifically we calculate the exponent a defined by means of the equation 

da^ da^ 

-ir = A ir< (10) 

where the superscripts refer to large nuclei (A) or deuterium (d) targets 

(3) 

and x 2 = x 2 . At this point it is instructive to observe that gluon fusion 
dominates the charmonium production process and is very important for 
bottomonium production as well. Therefore, the ratio of cross sections to a 
large extent reflects the ratio of gluon distributions, Rq^ . It is not hard to 
see that with the given definition of gluon distributions Rq monotonically 
increases with x 2 . As shown in Fig. 3, however, the data of Ref. |§ contradict 
this theoretical prediction. At x 2 ~ 0.15 there is a wide "bump" and the 
ratio starts decreasing at higher x 2 . The reason that this behavior becomes 
more apparent in the case of T production is the fact that the T is much 
more massive than the J/^f. For given center of mass energy, a particular 
quarkonium momentum, xf, probes a larger x 2 value in the T case. The 
gluon distributions being the most relevant and the least known among all 
the partons should be the first candidates for improvement. 



5.1 Improved gluon distributions 

At this point we turn to Information Theory. The total momentum fractions 
carried by the partons in an Nq cluster, 

z$= fdxFtfix), (11) 

<j 

where a designates the type of parton and F$ is the momentum distribution 
of a, must always add up to unity. Consequently the functions F$ satisfy 
the condition required in order to define an information entropy for them, 

S n = -£ C dxF { N a \x)\nF^ a \x). (12) 
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We shall maintain the quark distributions as defined in the previous sec- 
tion and modify the gluon distributions under the constraint that the total 
momentum fraction carried by gluons in each type of cluster is fixed and 
equal to that of the unmodified distributions, i.e., = 0.57 for nucleons 
and z^ = 0.60 for 6q clusters. 

The trend of the T data suggests that the simplest possible modification 
to the gluon distributions is an alteration of the linear term in x. We will, 
then, define corrected momentum distributions for the gluons in nucleons 
and 6q clusters as 

G N (x) = Cjv(1 - x) CN + C' N x. (13) 

For each N we now have two unknown parameters Cn and C' N . Momentum 
conservation, i.e., that the fractions z$ are constant, imposes the condition 

C' N = 2*<?> - ^L, (14) 
cn + 1 

where the exponents cn are kept fixed, C3 = 6 and cq = 10. Then we evaluate 
Cn from the requirement 

§^°- < 15 > 

The maximization procedure yields the following numbers: (C3, C' 3 ) = 
(1.163, 0.812) and (C 6 , C' % ) = (1.411, 0.987). We can compare these numbers 
with the uncorrected ones: (C 3 , C' 3 ) uncorr = (4.130, 0.0) and (C 6 , C' 6 ) uncorr = 
(6.624, 0.0). We have thus changed the shape of the function without affecting 
its integral, the total gluon momentum. The behavior of the new functions 
differs from that of the old ones mostly in the large x region for which on 
the other hand we have little experimental data. We note that the ocean 
distributions could not accommodate such type of alteration because deeply 
inelastic scattering imposes a constraint on the ratio of neutron to proton 
structure functions, F% \(x) / (x) — >l/4asx— > 1. 

Using the corrected distributions we can recalculate the exponent a and 
compare the results with the data. With a final state absorption cross section 
of 3 mbarn/nucleon and the new distributions the agreement with the data 
is considerably improved as shown in Fig. 3 in which the upper curve cor- 
responds to the uncorrected model, the lower curve to the corrected model 
using only gluon fusion contributions and the middle curve to the corrected 
model including gluon fusion and quark annihilation. It must be pointed 
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out that this is not the only model that gives reasonable description of the 



relative J/\I/ to T suppression data in p-A collisions. The authors of Ref. [30 
attribute the smaller suppression of T to the Q 2 ~ m 2 evolution of the distri- 
bution functions, where m is the mass of the resonance. Indeed, as discussed 



in Ref. [26 the evolution of the ocean (and consequently the gluon) distri- 
butions leads to smaller shadowing as Q 2 decreases. Our model neglects the 
Q 2 evolution relying on the observation [^(| that the masses of the quarko- 
nium resonances are already in the scaling region and attributes the reduced 
T suppression to its smaller absorption cross section. It is conceivable that 
both effects may, in fact, contribute to this phenomenon. We would mostly 
like to emphasize that it is the shape of the suppression curve that needed 
to be improved to account for the T data. 

It is worthwhile noting that due to the fact that T production probes a 
different kinematic regime from J / ip the correction on the gluon distributions 
has only a small effect on the J/\l/ suppression curves. The agreement with 
the charmonium data is still good as it can be observed in Fig. 4, although it 
slightly deteriorates at small xf (large X2). In addition the corrected model 
does not exhibit antishadowing of J/\l/ production for negative ij? a feature 
that would be in contradiction with the data PH . Data on bottomonium and 
charmonium give complementary information on the gluon distributions in 
nuclei. In addition the T is a really good probe of the initial state in which 
it is produced due to low absorption cross section. 



5.2 Predictions for RHIC 

We can use this model to make predictions for the J/^> and T suppression 
at the Relativistic Heavy Ion Collider (RHIC) with ^/s = 200 GeV/nucleon. 
The calculation proceeds along the same lines as for p-A collisions but now 
there is the additional possibility of 6q-6q cluster collisions. The nuclear 
effect is, thus, more pronounced. The cross section for collisions of a nucleus 
A with a nucleus B is 

d -^P- = KFd CZ £ E W *« + H ^ *n] . (16) 

dX F J4m c (b)2 1=3,6.7=3,6 

where = x^x^ /[m 2 (x^ )], i and j represent the type of colliding 
cluster, i/i is the nucleon-nucleon CM energy and H^JK are functions of 
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parton distributions appropriate for A-B collisions. 

In Fig. 5 we show theoretical results obtained with the uncorrected (curves 
marked by (d)) [^BJ and the corrected (curves marked by (c)) models. The 
curves exhibiting less suppression are for the T. The corrected model leads to 
larger suppression. The large Xp behavior is now much more distinct. It can 
be understood if we realize that in symmetric heavy ion collisions as xp — ► 1 
the large x region for the positively moving nucleus is probed. Therefore, in 
the uncorrected model the ratio increases with xp reflecting the increase in 
the gluon distribution ratio and exhibits antishadowing while in the corrected 
one it flattens out and remains below unity. 

6 Conclusions 

We have applied Information Theory to improve the gluon momentum dis- 
tribution functions in nuclei, including the "EMC effect" . The main idea is 
that by defining an information entropy, S, for those functions whose total 
integral is unity we can evaluate their parameters by maximizing S with re- 
spect to them. In other words we assume that the best choice of parameters 
is the one that is consistent with maximal ignorance under the constraint 
of momentum conservation. A quark-cluster model for the "EMC effect" 
has been used to establish good agreement with the data on J/ty suppres- 
sion in p-A collisions but proved inadequate to describe T suppression. Then 
Information Theory provided us with a tool to improve the model with signif- 
icant success. The gluon distributions have been corrected for their behavior 
at large x and an overall agreement with the differential cross sections for 
quarkonium suppression was achieved. 

An interesting aspect of this method is that it does not rely on any specific 
microscopic theory which in turn should be investigated in detail but uses 
only very general notions. The solution that is consistent with the assump- 
tion of maximal ignorance, quantified by the information entropy, seems to 
be an optimal one. This method can be used to improve the parameters of 
more detailed and realistic parton distributions in nucleons and nuclei under 
constraints imposed by experimental data. 

The author would like to thank R. Vogt, S. Gavin and A. Plastino for 
useful discussions. 
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Figure 1: The ratio of the experimental to the theoretical cross section for 
jj-ty production in collisions of protons with Cu and Be nuclei. The theo- 
retical cross section includes cluster contributions and final state absorption. 
The experimental cross section is from Ref. J7| with statistical and systematic 
errors added in quadrature. The uncertainty in the theoretical cross section 
is not included. 
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Figure 2: C (a), Ca (b), Fe (c), and W (d) to deuterium J/ty production 
ratio versus xp. The 6q probabilities are evaluated using / = k In A with 
0.040-0.052 for deuterium. The 6q ocean and gluon exponents are 10, 9, 
respectively for the solid lines and 12, 11 for the dotted ones. Panel (d) also 
shows the results for no absorption (short and long dash curves) and with 
absorption only (dot-dash line). The data are from Ref. || at */s = 38.7 
GeV. The errors are statistical. 
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Figure 3: Exponent of the W to deuterium T production ratio versus %2 ■ 
The chosen 6q probabilities are shown in the plot. The data points are from 
Ref. |J. The solid squares are for the 15 state and the open ones for 2S + 3S. 
The upper curve corresponds to the uncorrected model. The lower curve 
corresponds to the corrected model using only gluon fusion contributions. 
The middle curve corresponds to the full corrected model including quark 
annihilation. 
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Figure 4: The ratio of J/^f production cross section per nucleon on Cu to 
Be targets with a proton beam at \/s = 38.72 GeV. The solid curve is for 
the uncorrected gluon distributions and the dash curve for the modified ones. 
The exponents used are (c 3 , c 6 ) = (6, 10) for the gluons and (a 3 , a 6 ) = (9, 11) 
for the ocean. The effective 6q cluster probabilities are 0.30 for Cu and 0.16 
for Be. The data points are from Ref. 0. 
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Figure 5: The ratios of J/ty and T production Au-Au cross sections per 
nucleon over deuterium ones at RHIC energies. The curves marked by (d) 
are for the uncorrected model and those marked by (c) for the corrected one. 
The effective 6q probabilities are 0.40 for Au and 0.05 for deuterium. The 
used exponents are (c 3 ,c 6 ) = (6, 10) for the gluons and (a 3 ,a 6 ) = (9, 11) for 
the ocean. The absorption cross sections are 3.5 mbarn/nucleon for the J/^f 
and 3.0 mbarn/nucleon for the T. 
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